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The bound state wave functions for a wide class of exactly solvable potentials are found utilizing 
the quantum Hamilton-Jacobi formalism of Leacock and Padgett. It is shown that, exploiting the 
singularity structure of the quantum momentum function, until now used only for obtaining the 
bound state energies, one can straightforwardly find both the eigenvalues and the corresponding 
eigenfunctions. After demonstrating the working of this approach through a few solvable examples, 
we consider Hamiltonians, which exhibit broken and unbroken phases of supersymmetry. The nat- 
ural emergence of the eigenspectra and the wave functions, in both unbroken and the algebraically 
non-trivial broken phase, demonstrates the utility of this formalism. 

I Introduction 

Among st vari ous approaches to quantum mechanical problems JlrlllH the Hamilton-Jacobi 
formalismp^3 has been shown to be ideally suited for determining the bound state energies 
of exactly solvable problems. For a lucid review of various formulations of qua ntum mechanics and 
their relative merits and demerits, the reader is referred to Styer et. aiM^l In classical mechan- 
ics, the action-angle variables are the most appropriate ones for the description of the dynamical 
qq ' systems exhibiting periodic motional The corresponding demonstration, for t he qua ntum me- 
\Q . chanical bound state problems, is of recent origin. In 1984, Leacock and PadgetlPEl defined an 
exact quantum action variable through a quantum momentum function (QMF), whose quantization 
condition 

^ 1 i r 

^fpdx = nh, (1) 

' led to the bound state energies. The fact that only the residues of the QMF are required for 
Q^. finding the eigenvalues , and t he same can be evaluated rather easily from the quantum Hamilton- 
' Jacobi (QHJ) equatio n^ explains the economy of the QHJ formalism in finding the bound 
state spectra. In the limit h — ► the QMF goes over to the classical momentum, ^2m(E — V) 
and one can derive the well-known WKB quantization r ule. Th e demonstration of the connection 
between the supersymmetric WKB approximation schem ^^^ \ known to be exact for the so called 
shape-invariant potentials and the exact Hamilton-Jacobi qua ntization condition further clarified 
the working of various semi-classical quantization method*^. 

QHJ formalis m is fo rmulated in the complex plane, making it an idea l to ol for studying quasi- 
exactly solvabk P^^ and recently discovered PT-symmetric potentials^!, since these are also 
naturally studied in the same domain. It is well suited for studying chaotic systemJ23, where 
the periodic orbits are known to play a significant role. In this connection it is worth mentioning 
a rec ent demonstration of the close relationship between the Kolmogorov-Sinai entropy and the 
actional, both in classical and quantum domains. Before proceeding to the analysis of the above 
mentioned non-trivial quantal problems, it is necessary to develop a precise understanding of the 
procedure for finding the wave functions from the quantum Hamilton-Jacobi equation. Further- 
more, in light of the close similarity between the QHJ and the WKB type quantization conditions 
and the well - known difficulties one encounters in extracting the wave functions in the semi - 
classical approaches, it is of deep interest to find the same in the QHJ scheme. 

The present work is devoted to the study of bound state wave functions through the quantum 
Hamilton-Jacobi formalism. It is shown that the singularity structure of the quantum momentum 
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function, which has been so far used only for finding the bound state energies, can straightforwardly 
yield both the eigenvalues and the corresponding eigcnfunctions. In the earlier works of, two of the 
present authors and those of Leacock and Padgett, the singularity structure of the QMF and the 
boundary conditions it needs to satisfy have been carefully studied. Since the same potentials are 
being analyzed here, we use the singularity structure of the QMF found earlier to determine the 
eigenstates. We will al so not deal with the details of the QHJ formalism and refer the interested 
reader to earlier paperfJlUl^l. 

We first demonstrate the working of this approach through a number of exactly solvable po- 
tentials, which have been studied extensively in the literature thro ugh a number of approaches, 
ranging from factorizatio n m ethooP^Bl path integral techniqueJ^El, group theoretical method^ to 
SUSY operator methocPJEH. 

We then analyze Hamiltonians which exhibit broken and unbroken 
phases of supersymmetry. The natural emergence of the eigenspectra and the wave functions, 
in both the phases of supersymmetry, demonstrates the utility of this formalism. It should be 
mentioned that Hamiltonians exhibiting broken phase of supersymmetry are not straightforwardly 
amenable to an algebraic treatment. 

The paper is organized as follows. We first give a brief description of the QHJ formalism 
using the well-studied harmonic oscillator problem as a prototype. We explicate the procedure 
to obtain the bound state wave functions and its connection wit h th e singularity structure of 
the QMF. This is an example of the confluent hypergeometric clasd^l. We then proceed to the 
study of Rosen-Morse and the Scarf - I potentials, belonging to the hypergeometric class. It is 
interesting to note that, Scarf-I potential, although similar to the Rosen-Morse, exhibits both 
broken and unbroken phases of supersymmetry, for different ranges 

of potential parameter^ 

The role of the proper implementation of the boundary conditions, satisfied by the QMF in 
Hamilton- Jacobi formalism, is clearly brought out in this section. It is shown that, imposing the 
boundary conditions correctly leads automatically to the appropriate bound state wave functions, 
for the broken and the unbroken phases of supersymmetry (SUSY). 

II Harmonic Oscillator in the QHJ formalism 

We begin this section with a brief summary of the QHJ formalism for one dimensional 
potentials. As pointed out by Leacock and Padgett, the QHJ equation, 

p 2 -ih^ = 2m(E-V(x)), (2) 
ax 

is a special case having the well-known Riccati equation, 

A(x) P 2 +B(x) P + C{x) + i J = 0. (3) 
@ can be connected to the Schrodinger equation for a given potential V(x) as 

h 2 (Pip(x) 



(/)J +V(x)iP(x)=E 1 p(x), (4) 
provided the QMF p is related to the wave function ip(x) by 

p = — ih— \nip(x). (5) 
ax 

The above relation can be inverted to yield tp(x) — exp (i J pdx). We first analytically continue p 
to the complex plane by taking x as a complex variable. In what follows, the QMF plays a central 
role in finding the bound state wave functions. The knowledge of the singularities of the QMF 
is sufficient to determine its form, which conseq uen tly lead to the wave functions. The solutions 
of Riccati equation has two types of singularitieJ^HI ( the fixed and the moving singularities. The 
fixed singular points of p(x) reflect the singularities of the potential and are energy independent. 
Besides these, the QMF has other singular points called the moving singularities of the type, 

P~7Z ( 6 ) 
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which cancel amongst the terms in the left hand side of (01. It can be seen clearly from (|2J|, that 
the residue is given by — ih. These poles, unlike the fixed singularities, are dependent on energy 
and the initial conditions. It is known that, for the solutions of the Riccati equation, a moving 
singularity can only be a polcP^2EI3 it is a well- known fact that, the wave function of the n th 
excited state will have n zeros on the real line in the classical region, which implies that the QMF 
will have n poles. In all the solvable models studied within the QHJ approach, no other moving 
singularity was found excep t the n moving poles. This observation has been very useful in solving 
for the bound state energie^^EI i n a u the potentials analyzed here, it has been assumed that, 
the QMF has finite number of singularities in the complex plane and that the point at infinity is 
an isolated singular point. In the present study we find that, after a suitable change of variable 
followed by transformations on p to bring the QHJ equation to the form of (J2J, the QMF is a 
rational function for the potentials studied here. The knowledge of the location of poles and the 
corresponding residues are sufficient to fix the form and to obtain the bound state wave functions. 

The QHJ equation can be used to obtain the residues at different poles. Owing to the fact 
that the Riccati equation is quadratic, one gets two answers for the residues. Hence, one needs a 
boundary condition to select the correct value. The original boundary condition 

limp— >Pd, (7) 

which was given by Leacock and Padgetl P^ I ^ I can be cast in several equivalent forms, which will 
be given as and when they are used. We will use one of the convenient forms to identify the correct 
residue. 

For the purpose of illustration, we start with the harmonic oscillator V(x) — l/2(kx 2 ). The 
Schrodinger equation, with £ = ax, A = 2E/haj c ,uj c — ^Jk/m and mk = % 2 a A , 

^ + (A-a^)=0, (8) 
has the corresponding QHJ equation in the form 

✓ -<!-*-<•■ m 

From the above equation it appears reasonable to assume that for large £, p goes as ±i£. To fix the 
right behavior of p at infinity one can use the boundary condition Q or instead use the relation 
between p and tp 

ip(x) = exp(i J pdx). (10) 

Demanding the square integrability of the wave function, the correct behavior for large £ is seen 
to be 

lim p i£. (11) 

The QMF p has n moving poles on the real line with residue — i. There are no fixed singularities 
in the finite domain in the complex plane, since the potential is not singular. It is p ossible to show 
that there are no other moving poles of the QMF in the finite complex plandHH^l thus making p 
analytic everywhere except at the n moving poles. Therefore we can write p as 

n 

p(0 = £ttV + ^ + Q(0, ( 12 ) 

k=l 5 ^ k 

where the summation term describes the sum of all the principal parts of the individual Laurent 
expansion of p around each moving pole. Q(£) is an entire function bounded for large £, with i£ 
describing the behavior of p at oo. Thus from Liouville's theorem Q(£) is a constant C. Hence the 
expression (|12|l can be cast in a useful form 

p(£) = -i^r+it + C (13) 
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where P is an n th degree polynomial IIfc=i(£ — £fe) an d ^ s seen to satisfy 

P'^ - 2CP r ; - 2CiP' n - (2i£C + C 2 - A + 1)P„ = 0. (14) 

Equating the coeffecients of £ n+1 and £" to zero gives 

C = 0, A = 2n+ 1 (15) 

and hence P„ = (n + ^)Tiu! c . With these values 114f> becomes the Hermite equation and P n 
coincides with Hermite polynomial H n (^) apart from an overall constant and one obtains the 
known eigenfunctions 



V>(£) = CX P (i J P d £ 



= exp(-|-+ln(ff n (0) (16) 

One can solve the radial part of hydrogen atom in the same way as the above problem except 
that one encounters a fixed pole at the origin along with the n moving poles. Thus one needs to 
take this pole into account, when one writes the behaviour of the QMF of the hydrogen atom in 
the entire complex plane similar to (|15J) . For more details one is referred to which is an earlier 
version of this paper. The bound state wave functions for other solvable potentials, belonging to 
confluent hypergeometric case, e.g., Morse and the radial oscillator potentials can be obtained in 
the same way. In some cases like Morse potential, one needs a suitable change of variable, before 
solving the QHJ equation. 

Ill Rosen-Morse Potential 

This section is devoted to a study of solvable potentials belonging to the hypergeometric 
class. Note that, in general one does not come across potentials which are meromorphic in nature. 
As mentioned earlier to make the QMF a rational function one needs to do a suitable change of 
variable y = f(x). For a general potential the QHJ equation 

p 2 - ip' -{E- V(x)) = (17) 

after a change of variable y = f(x) and introduction of q by q(y) = ip(x(y)), becomes 

q 2 (y) + F(y)^±+E-V(y) = 0. (18) 
ay 

Here F(y) equals the derivative expressed as a function of y and V(y) = V(x(y)). We see 

that l|18f) does not have the convenient form of ©. To bring (|18|l to the QHJ form we introduce 
X by means of the following equations. 

q = F{y)<t>, X =<t>+~{^F{y)), (19) 
2 dy 

and obtain 

X + dy + (F(y)r 2 { F(y) J + a{ F{y) ) ( 

One can observe from H2UI) . that the residues at the n moving poles of x will be unity. Therefore, 
we can make use of l|20|) . instead of (|1 7|) . for any general potential and proceed in the same way 
as for the harmonic oscillator. 

For the Rosen-Morse potential, where the expression for the potential is 



V(x) = A 2 - A(A + a)sech 2 ax 



(21) 
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with A > 0, we first change the variable to y = tanh(ax). Using l|19(l . one obtains the equation for 
X as follows : 



o J5 - A 2 + a 2 ALA + a) 



X + + a 2 (l-y 2 ) 2 + a 2 (l-y 2 ) °' ^ 

For the n th excited state, besides having n moving poles with residue equal to one, QMF will have 
fixed poles at y — ±1. We assume that there arc no other singular points in the finite complex - y 
plane. Thus proceeding in the same way as in the harmonic oscillator, one can write x(y) as the 
sum of the principal parts of the individual Laurent expansions around different singular points, 
plus an analytic part which will turn out to be a constant C. Thus 

*>^ + ^T + ? + C; (23) 

here b\ and b\ are the residues at y = ±1 respectively. One obtains the residues at these poles by 
doing a Laurent expansion of x{y) around each pole and substituting it in the QHJ equation for 
x(y)- This leads to the values 



1 \J A 2 — E , 1 VA 2 -E 

6i = -±- and b =-±- . (24 

2 2a 1 2 2a y ' 

The ambiguity in the signs of the square roots can be eliminated by applying the boundary con- 
dition 0). Since we are dealing with the SUSY potential, we can make use of the boundary 
condition, 

lunq(y) = -W(y), (25) 

hj — >u 

instead of the boundary condition Q proposed by Leacock and Padgett, with W — Ata,nh(ax), 
the superpotentiaM Thus one obtains the correct residues as follows 



1 VA 2 -E ,. 1 VA 2 -E 
6l =2 + ^^ 6i =2 + ^^- (26) 

Substituting x{v) with the right values of b\ and b[ in (|22l) . and for large y comparing the coefficients 
of 1/y and the coefficients of 1/y 2 , one gets 

C = 0, E = A 2 - (A - naf (27) 

respectively and is left with the differential equation 



2P'y a + VA 2 -E \ {a + V A 2 - E) 2 E - A 2 + a 2 A(A + a) 
~P + ~~P~ I a(y 2 - 1) J + 2a 2 (y 2 - 1) 2a 2 (y 2 - 1) ~ a 2 (y 2 -l) ~ ' ( ' 

Substitution of the energy eigenvalue in the above equation gives the Jacobi differential equation 
for P 

(1 - y 2 )P" -y{n + v+ 2)P' + n(fx + v + n + l)P = Q. (29) 

Thus P is a Jacobi polynomial represented by P£ ,l> {y), where fj, = v — — — n. The bound state 
wave functions can be obtained by doing the change of variable and using , in (fLljl gives 

iP(y) = exp (J x (y)dy + , (30) 

where x{y) i s 



. a + VA 2 - E a + ^A 2 -E P> 

xiy) = Mv-i) Mv + i) + P' ( } 
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Thus the bound state wave function from l|30|) turns out to be 

tfn(y) = JV(i-y)*(i + v)*J*r(y) ( 32 ) 

which matches with the known results "3 for A, real and positive. 
IV Scarf - I (trigonometric) 

There exist 

potentialfPH2Hl 

for which SUSY can be exact or broken, depend ing on the range 
of potential parameters. As mentioned earlier in the introduction, Bhalla et aly^ have studied 
the eigenspectra of Scarf - I potential. In this section, we study this potential and obtain the 
expressions for the wave functions. The expression for the potential is 

V_ (x) = -A 2 + (A 2 + B 2 - Aah) sec 2 ax - B(2A - ah) tan ax sec ax. (33) 

For this potential, SUSY is known to be exact in the parameter range 

(A — B) > 0, (A + B) > 0, (34) 

and broken in the range 

(A - B) >0, (A + B) <0. (35) 



Since broken SUSY is one of the possibilities, unlike the previous cases, we cannot assume SUSY to 
be exact and hence cannot use the superpotential to find the correct residue at the fixed pole s. We 
therefore, fall back upon the boundary condition originally proposed by Leacock and Padgetl 



mm 



limp^pd. (36) 

It turns out that imposing the boundary condition consistently for different ranges of the param- 
eters A and B, automatically gives rise to the correct wave functions and energy eigenvalues for 
both broken and unbroken phases of SUSY. Retaining ft with (2m = 1), replacing p with —iq in 
the QHJ equation, one gets, 

q 2 + h^- + E + A 2 - ( A 2 + B 2 - Aah) sec 2 ax + B(2A - ah) tan ax sec ax = 0. (37) 
ax 

Making a change of variable 

y = sin ax, (38) 

and switching to \ in place of q as explained previously, we get 

2 n dx V 2 h 2 E + A 2 a 2 h 2 - 2(A 2 + B 2 - Aah) B(2A - ah)y 

X dy 4(1 -y 2 ) 2 a 2 {l-y 2 ) 2a 2 (\-y 2 ) 2 a 2 {\-y 2 ) 2 [ ' 

Proceeding in the same way as before, we assume that there are no singular points in the finite 
complex plane, other than the fixed poles at y = ±1 and the n moving poles on the real line. 
Therefore, we write x(u) as 

x(y) = -^T + -XT + ^ + c, (40) 

y-l y+1 P 

where b\ and b[ are the residues at y = ±1 respectively. The values of the residues bi and b\ are 
obtained as in the previous case. The values of b\ at y = 1 are 

(A-B) h , (A - B) 3ft 
The two values of the residue at y = — 1 are 
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Different pairs of b\ and b\ give physically acceptable solutions in the two phases of SUSY. We 
proceed to obtain these solutions by taking the correct choice of the residues using (|3(jf) . the 
details of which are similar to those in Bhalla ei.aZp^l 1] The boundary condition for this range 
of values of A and B implies, b\ > and b[ > in the limit U — > 0. 

Case:l A + B > 0,A- B > 

For this particular range of parameters the right choice of the residues is 

Thus x m l|^fJ|) becomes 

. . (A-B h\ 1 fA + B h\ 1 P> 

x{y)= [— + ±) — + {— + VVTT + T (44) 

substitution of which in 139|) and putting h = 1 gives the energy eigenvalue 

E n = {A + naf - A 2 . (45) 

The equation for P assumes the form 

(1 - y 2 )P" + + l)y)P' + n(n + fi + v)P = 0, (46) 

where \x — A ~ B and v = A ~^ B . The above differential equation is of the form of the Jacobi 

differential equation and hence P coincides with P„ 1 ^ 2,l/ 1 ^ 2 . The bound state wave function 
when SUSY is exact, is 

My) = N(l - y)*Q. + y)*i^~*'" _4 (|/). (47) 
These wave functions correspond to the solutions of the unbroken phase of SUSY.® 
Case:2 A- B > 0,A+ B <0 

For this range of the parameters, the right choice of the residues is 

(A-B) h , (A + B) 3H 

Substitution of \ with these values of b\ and b' x in (|39(l and putting h — 1 gives the energy 
eigenvalue as 

1, x2 



E= IB- („+-)«) - A , (49) 

and the differential equation for P turns out to be 

(1 - y 2 )P" + (0 - ii + 1) - y{v + n + 2))P' + n(n + 1 + n + v)P = (50) 
and the bound state wave function is found to be 

My) - iV(i - y)^ (i + y)*-*Pn +hfl ~Hy) (51) 

where These bound state wave functions correspond to the broken phase of 

SUSY of the Scarf-I potential. The Scarf-I potential can be related to the Poschl-Teller-I potential 
by redefinition of the potential parameters as 

ax ^ ax + A = S + (3 and B = S — (3. (52) 
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With this redefinition of the parameters and then putting a = 1, we see that the bound state 
wave functions in (|5 1|) matches with the bound state wave functions of the broken SUSY phase of 
the P6schl-Teller-I potential given i Thus we see that the QHJ formalism in one dimension 
gives the correct expressions of the bound state wave functions, when there are different phases 
of SUSY. It may be remarked here that, in the range A — B < 0, A + B < 0, SUSY is exact but 
the roles of H_ and H + are interchanged. In the range A — B < 0,A + B > SUSY is again 
broken. The methods outlined above, has been applied to many other potentials like the Morse, 
Eckart, P6schl-Tellcr(hyperbolic) and Scarf-II (hyperbolic) and the results agree with those given 
inP 

V Conclusions 

In conclusion, the singularity structure of the quantum momentum function, in the complex 
domain, can be effectively utilized for obtaining both eigenfunctions and eigenvalues in an efficient 
manner. Besides illustrating the procedure to find the wave functions, through the analysis of 
the exactly solvable problems, we have clarified the role of an appropriate implementation of 
the boundary condition, necessary to solve the quantum Hamilton- Jacobi equation. The same 
played an important role in yielding the eigenvalues and eigenfunctions of potentials exhibiting 
broken and unbroken phases of SUSY. As has been pointed out earlier, a close connection 
between Kolmogorov-Sinai entropy and the action has been obtained recently, making use of a 
Hamilton- Jacobi approach through Bohm's formulation The QHJ formalism of Leacock and 
Padgett and the above are not identical. Hence, a deeper exploration of the above connection, 
as well as the extension of the Hamilton- Jacobi formalism to intrinsically non-separable systems, 
like Hamiltonians whose classical dynamics reveal chaotic behavior, need thorough investigation. 
Although, we have dealt with solvable systems here, one needs to extend this approach, to 
non-solvable cases, like various anharmonic oscillators. We would like to mention that, during 
the course of our study, we have analyzed scattering states of the Rosen-Morse potential. It was 
found, by assuming that QMF has finite number of moving singularities, one can obtain the 
scattering state wave functions, provided the potential parameter A is equal to an integer n, a 
special case when the potential is reflectionless. The case when A^n needs further study, as also 
a more careful analysis of the general scattering problem. The continuous energy functions, of 
interest for scattering, are currently under study. 
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